Molecular signatures in the structure factor of an interacting Fermi gas. 
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The static and dynamic structure factors of an interacting Fermi gas along the BCS-BEC crossover 
are calculated at momentum transfer ftk higher than the Fermi momentum. The spin structure factor 
is found to be very sensitive to the correlations associated with the formation of molecules. On the 
BEC side of the crossover, even close to unitarity, clear evidence is found for a molecular excitation 
at h 2 k 2 /4m, where m is the atomic mass. Both quantum Monte Carlo and dynamic mean-field 
results are presented. 
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The possibility of producing weakly bound molecules 
in interacting ultracold Fermi gases, raises very interest- 
ing challenges. These molecules are formed near a Fesh- 
bach resonance for positive values of the s-wave scattering 
length, have bosonic character and exhibit Bosc-Einstein 
condensation at low temperature 0, Q . They have a re- 
markably long life-time as a consequence of the fermionic 
nature of the constituents which quenches the decay 
rate associated with three-body recombinations. Several 
properties of this new state of matter have been already 
investigated experimentally in harmonically trapped con- 
figurations. These include the molecular binding en- 
ergy 0. the release energy , the size of the molecular 
cloud 0] , the frequency of the collective oscillations [j| , 
the pairing energy Q , the vortical configurations 0] and 
the thermodynamic behavior p|. 

In this paper we investigate another feature of these 
new many-body configurations, directly related to the 
molecular nature of the constituents: the behavior of 
the static and dynamic structure factor at relatively high 
momentum transfer. Experimentally the structure factor 
can be measured with two-photon Bragg scattering where 
two slightly detuned laser beams are impinged upon the 
trapped gas. The difference in the wave vectors of the 
beams defines the momentum transfer ftk, while the fre- 
quency difference defines the energy transfer fku. The 
atoms exposed to these beams can undergo a stimulated 
light scattering event by absorbing a photon from one of 
the beams and emitting a photon into the other. This 
technique, which has been already successfully applied 
to Bose-Einstein condensates @, provides direct access 
to the imaginary part of the dynamic response function 
and hence, via the fluctuation-dissipation theorem, to the 
dynamic structure factor. At high momentum transfer 
the response is characterized by a quasi-elastic peak at 
lu = fik 2 j2M, where M is the mass of the elementary con- 
stituents of the system. The position of the peak is conse- 
quently expected to depend on whether photons scatter 
from free atoms (M — m) or molecules (M = 2m). For 
positive values of the scattering length both scenarios 
are possible and their occurrence depends on the actual 
value of the momentum transfer. If k is much larger 
than the inverse of the molecular size, photons mainly 



scatter from atoms and the quasi-elastic peak takes place 
at Hk 2 /2m. In the opposite case, photons scatter from 
molecules and the excitation strength is concentrated at 
hk 2 1 Am. The two regimes are associated with different 
velocities of the scattered particles given, respectively, by 
hk/m and hk/2m. 

Let us suppose that the Fermi gas consists of an equal 
number N/2 of atoms in two hypcrfme states (here- 
after called spin-up and spin-down). Using Sfj (k, u>) — 
Sn(k,uj) and S^(k,u>) = S^(k,u)), we can write the 
T = dynamic structure factor in the form 



S(k,L0) = 2(S n (k,L0)+S tl (k,L0)) , 



(1) 



with 



S a „<(k,u))=^2< 0\p a (k)\n >< n\pl,(k)\Q > 5{hu-E n0 ) 

n 

(2) 

where p a {k) — J2i e %kZi are the spin-up (a =|) and 
spin-down (a =1) components of the Fourier transform of 
the atomic density operator, while \n > and E n o = E n — 
Eq are the eigenstates and eigenenergies of the many- 
body Hamiltonian H . 

The frequency integral of the dynamic structure factor 
defines the so-called static structure factor relative to the 
different spin components: 
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The total spin structure is then given by S(k) = 
N^h JdujS{k,uj) = S n {k) + S n (k). The static struc- 
ture factor is related to the two-body correlation func- 
tions through the relationships 



S TT (fc) 
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yielding S(k) = 1 + n J dr [g(r) — l]e~ lk ' r with g(r) — 
[#Tt( r ) + 5ti( r )]/2- In the above equations n is the total 
particle density fixing the Fermi wave vector according 
to kp = 3ir 2 n. The behavior of the structure factor S(k) 
at small momenta {k -C kp) is dominated by long-range 
correlations which give rise to a linear dependence in k. 
In a superfluid the slope is fixed by the sound velocity 
c, through the general law S(k) = k/2mc. In this paper 
we are however mainly interested in the behavior at large 
momentum transfer, typically such that k > kp. 

In the limit of very large k, the sum in Eq.Q is dom- 
inated by the autocorrelation term i = j with identical 
spins. This leads to S^(k) — » 1. On the other hand, 
since there is no autocorrelation with different spins, 
Sn(k) — > for very large k. Actually in the ideal Fermi 
gas the dynamic spin-up spin-down structure factor iden- 
tically vanishes (S^i(k,uj) = 0) for all values of k and w, 
reflecting the complete absence of correlations between 
particles of opposite spin |lOj| . This quantity is therefore 
particularly well suited for studying the effect of interac- 
tions. Let us first consider the case of small and positive 
scattering length kpa <C 1. This is the so-called Bose- 
Einstein condensation (BEC) regime, where we have a di- 
lute gas of weakly bound molecules, made of atoms with 
opposite spins, with normalized wavefunction <I>o( r ) for 
the relative motion. When we consider distances of the 
order of the molecule size a, we have naturally a strong 
correlation between opposite spin atoms belonging to the 
same molecule. In this case the sum is dominated by 
this contribution, which gives: 

S n (k) = J dv e^ r n mol {v) , (6) 

where n mo i(r) = |$o( r )| 2 is the probability to have the 
atoms separated by r. This holds only for k 3> fcp , oth- 
erwise one has to take into account also correlations be- 
tween atoms belonging to different molecules. In partic- 
ular one finds that, for kp <C k <C 1/a, S^i(k) ~ 1, so 
that S(k) = 2. This corresponds to the regime where the 
elementary constituents "seen" by the scattering probe 
are molecules and not atoms. 

When one moves away from this BEC regime towards 
the resonance, where kpa 3> 1, the many-body wave 
function cannot be simply written in terms of molecules 
anymore. In this very interesting regime the function 
S-\i(k) is smaller than unity, but can still significantly 
differ from zero, reflecting the crucial role played by in- 
teractions (see Fig. Pi . 

In Figs. 11131 we report quantum Monte Carlo (QMC) 
results of the static structure factors as a function of 
k/kp, calculated for different values of the relevant inter- 
action parameter kpa. The simulations are carried out 
following the method described in Ref. [ll|. The direct 
output of the calculation are the spin-dependent pair cor- 
relation functions, g-\-\(r) and gn(r). The corresponding 
structure factors are obtained through the Fourier trans- 
formation of Eq. JSJ. Notice that the low- momentum 
behavior of the structure factors can not be accessed in 
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FIG. 1: (color online). Spin-dependent static structure factor 
S]](k) and Sn(k) (inset) obtained using the QMC method, 
for different values of the interaction strength as indicated in 
the figure. The red dotted line in the inset corresponds to the 
Fourier transform of n mo ;(r) [Eq. 




FIG. 2: (color online). Total static structure factor S(k) for 
different values of the interaction strength. Solid lines corre- 
spond to QMC results and dashed lines to BCS results. The 
black line refers to the non-interacting gas. 

the QMC simulation due to the finite size of the simula- 
tion box. For our system (iV=66 particles and periodic 
boundary conditions) a reliable calculation of S aa i (k) is 
limited to k > 0.5 kp. 

The results for Sf [(k) for small positive values of kpa 
(see Fig. ^ are fully consistent with the transition from 
the molecular regime, where S-\i(k) ~ 1, to the atomic 
one where Sjj(fc) — > at large k. The transition is 
well described by Eq. JHJ which is reported in the in- 
set of Fig. n for the value kpa = 1/4. The molecular 
density profile n mo /(r) entering Eq. JfJ) has been cal- 
culated using the same two-body potential of range Rq 
employed in the QMC simulation jl2| . In the large-fc re- 
gion, l/Ro ^> k ^> max(kp, 1/a), the decay of Sn(k) is 
proportional to 1/fc as a consequence of the 1/r 2 short- 
range behavior of .0||(r). This behavior is fixed by two- 
body physics which dominates at short distances. The 
proportionality coefficient of the 1/fc law is given by 
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FIG. 3: (color online). Magnetic static structure factor Sm (k) 
for different values of the interaction strength. Solid lines 
correspond to QMC results and dashed lines to BCS results. 
The black line refers to the non-interacting gas. 

27T 2 lim r .^ [r -2 1 <E>o (r) | 2 ] in the deep BEC regime (r — » 
here means Rq <C r <C a) and by many-body effects closer 
to the resonance and in the BCS regime. For example, 
at unitarity (1/kpa = 0), we find Sn(k) ~ 0.85fcp/fc. 
For 1/kpa 1, the "atomic" regime, where Sn(k) ap- 
proaches zero, is reached only at extremely large values 
of k/kF reflecting the crucial role played by interactions 
on this quantity even at high k. Conversely is less 

sensitive to interactions unless one considers small values 
of k. This quantity approaches quite rapidly its large- A: 
limit St T (fc) -> 1 (see Fig. 0. 

The results for the total structure factor are shown 
in Fig. |21 In the BEC regime S(k) first increases and 
reaches a plateau where S(k) ~ 2, which is characteristic 
of the molecular regime. At higher k it decreases and 
eventually approaches the uncorrelated value S(k) = 1. 
In Fig. |3 we show the results for the magnetic structure 
factor Sjwr(fe) = Sjj(fc) — Sfi(k). Compared to the to- 
tal structure factor this quantity in the BEC regime is 
strongly quenched in the range ka <C 1 where it behaves 
like (ka) 2 . In fact, the leading contributions of S^(k) 
and S-\i(k), of order of 1, cancel each other in this regime. 

In Figs. 12131 we have also reported the results of the 
dynamic mean-field theory (see discussion below) which 
provides very reasonable predictions for the structure fac- 
tors as appears from the comparison with the microscopic 
findings of the QMC simulation. The deviations at large 
k in the BEC regime are mainly due to the details of the 
molecular wave function used in the QMC simulation at 
distances r < R . 

Let us now discuss the behavior of the dynamic struc- 
ture factor. Since QMC calculations for this quantity 
are not available, we will discuss here the predictions 
of the dynamic self-consistent mean-field approach based 
on the Bardeen-Cooper-Schrieffer (BCS) theory, which 
has proven to be quite accurate in the evaluation of the 
static structure factors (see Figs. [3 El ■ T° evaluate the 
dynamic structure factor in the self-consistent BCS the- 
ory a convenient procedure is given by the calculation of 



FIG. 4: (color online). Dynamic structure factor, from the 
BCS approach, as a function of uj for fixed momentum transfer 
k = 3&f and for various interaction strengths. Blue solid lines 
refer to the total structure factor S(k,u) and red dashed lines 
to the magnetic one SW (fc, uj). The black dotted line in panel 
d) corresponds to the non-interacting gas. In panel a) and b) 
is also indicated the weight of the sharp peak to S(k). The 
units of the dynamic structure factor are N/(2tF). 

the linear response function. This is obtained by writing 
a kinetic equation which generalizes the well known Lan- 
dau equation for Fermi liquids to the superfluid state. 
This equation includes as driving terms the off-diagonal 
self-consistent field, which ensures that the resulting the- 
ory is gauge invariant, in contrast to the original ele- 
mentary BCS theory. In particular this approach leads 
naturally to the appearance of the gapless Bogoliubov- 
Anderson mode on the BCS side of the crossover [l3j . 
which goes continuously into the Bogoliubov mode in 
the opposite BEC regime On the BCS side this 

collective mode exists only at low wavevector k. When 
k is increased it merges at some point into the contin- 
uum of single-particle excitations. By approaching the 
Feshbach resonance from the BCS side, the threshold of 
single-particle excitations is pushed towards higher ener- 
gies, which makes the merging occur for higher k. Rather 
soon, on the other side of the resonance, i.e. for small 
positive 1/kpa, the mode dispersion relation never meets 
the continuum anymore, and it evolves to the standard 
Bogoliubov mode of the molecular BEC regime, the con- 
tinuum of the single-particle excitations being located at 
higher energies 15]. At the high momenta considered in 
this work the Bogoliubov mode takes the free molecule 
dispersion ui(k) = hk 2 /Am. This result is fully consis- 
tent with the behavior of the static structure factor S(k) 
discussed above. In fact using the f-sum rule 



dio u)S(k, uj) = N 



h 2 k 2 

2m 



(7) 



and assuming, following Feynman, that the integral is 
exhausted by a single excitation at uj(k), one finds 
hui(k)S(k) — h 2 k 2 /2m, which yields the dispersion 
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oj(k) = Hk 2 /4m in the molecular regime where S(k) = 2. 

The possibility of seeing these "free molecules" in the 
dynamic structure factor depends in a critical way on the 
value of the dimensionless combination ka. In Fig. 0] we 
show the behavior of the dynamic structure factor as a 
function of fkj/ep, where ep — h 2 k 2 F /2m is the Fermi en- 
ergy, for a fixed value of momentum transfer (k/kp = 3), 
so that the energy-weighted integral 0) of the dynamic 
structure factor takes the same value for all the configura- 
tions in the figure. Different values of kpa are considered. 
When kpa is positive and small the system is in the BEC 
regime. Fig.^Jt shows that in this case (1/fc^ a = 2) most 
of the strength is concentrated in the molecular peak lo- 
cated essentially at oj(k) — Kk 2 /Am (corresponding to 
frui/cp = 3 2 /2 = 4.5). The shift of the position of the 
peak with respect to the "free" molecule at Hk 2 /4m re- 
flects interaction effects between the molecules. To low- 
est order it is given by the mean-field result Scu = gn/2h, 
where g = 4irh 2 a m /2m is the molecular coupling con- 
stant with a m denoting the molecule-molecule scattering 
length (a m — 2a in the BCS approach, but a m = 0.6a in 
exact treatments 0). For l/k F a = 1 the strength in the 
continuum of atomic excitations is larger, but the molec- 
ular peak is still quite strong since it carries 89.5% of the 
total strength. At unitarity (Fig. the peak has just 
merged into the continuum of the single-particle excita- 
tions, whose spectrum turns out to be strongly affected 
by molecular- like correlations. Finally if one considers 
the BCS regime of negative k F a (Fig. 0Ji) one finds that 
rapidly the response is pratically indistinguishable from 
the one of the ideal Fermi gas and no residual effect of 
interactions is visible at such values of momentum trans- 
fer. 

In the same figure we have also shown the prediction 
for the magnetic dynamic structure factor SM(k,uj) — 
2[Sii(k, u) — S^i(k,uj)]. Compared to S(k,uj) this quan- 
tity is not affected by the presence of the collective mode. 
Furthermore its strength in the continuum turns out to 
be pushed up to higher frequencies. This is consistent 
with the fact that SM{k,u) obeys the same f-sum rule 



(0) as the total structure factor S(k,ui). 

The calculations presented in this work hold for a uni- 
form gas. At the high momentum transfers considered 
here, the effects of trapping, relevant for actual experi- 
ments, are expected to give rise to a small Doppler broad- 
ening in the molecular peak on the BEC side, similarly 
to what happens in Bose-Einstein condensates 0. As 
far as the static structure factor is concerned, the main 
effects due to the non uniform confinement can be taken 
into account using a local density approximation. 

As already anticipated the dynamic structure factor 
can be measured through inelastic photon scattering. In 
order to have access separately to the total and spin com- 
ponents of the structure factor various strategies can be 
pursued. A first possibility is to adjust the detuning of 
the laser light in such a way as to produce a different 
coupling with the two spin components, resulting in a re- 
sponse of the system that will also be sensitive to the spin 
structure factor. Another possibility is to work with po- 
larized laser light which can result in a different coupling 
in the two atomic spin components even if the detuning 
is large [Tt| . 

In conclusion we have shown that the structure fac- 
tor of an interacting Fermi gas at relatively high mo- 
mentum transfer exhibits interesting features associated 
with the molecular degrees of freedom of the system. In 
particular in the BEC regime and for values of k in the 
range kp < k < I/a, the total static structure factor ap- 
proaches the value 2, characteristic of a molecular regime, 
while the dynamic structure factor is characterized by the 
occurrence of a sharp peak at u>(k) = Hk 2 /4m. We have 
found that the up-down spin component of the structure 
factor exhibits a specific dependence on the momentum 
transfer, being strongly sensitive to the correlations asso- 
ciated with the formation of molecules. Bragg scattering 
experiments to explore the structure factors close to res- 
onance would be highly valuable. 
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